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Abstract 

Two-dimensional almost-Riemannian structures are generalized Riemannian structures on 
surfaces for which a local orthonormal frame is given by a Lie bracket generating pair of vector 
fields that can become coUinear. Generically, the singular set is an embedded one dimensional 
manifold and there are three type of points: Riemannian points where the two vector fields 
are linearly independent, Grushin points where the two vector fields are collinear but their Lie 
bracket is not and tangency points where the two vector fields and their Lie bracket are collinear 
and the missing direction is obtained with one more bracket. Generically tangency points are 
isolated. 

In this paper we study the Laplace-Beltrami operator on such a structure. In the case 
of a compact orientable surface without tangency points, we prove that the Laplace-Beltrami 
operator is essentially self-adjoint and has discrete spectrum. As a consequence a quantum 
particle in such a structure cannot cross the singular set and the heat cannot flow through the 
singularity. This is an interesting phenomenon since when approaching the singular set (i.e. 
where the vector flelds become collinear), all Riemannian quantities explode, but geodesies are 
still well defined and can cross the singular set without singularities. 

This phenomenon appears also in sub-Riemannian structure which are not equiregular i.e. 
in which the grow vector depends on the point. We show this fact by analyzing the Martinet 
case0 



^This research has been supported by the European Research Council, ERC StG 2009 GeCoMethods, contract 
number 239748, by the ANR Project GCM, program Blanche, project number NT09-504490 and by the DIGITEO 
project CONGEO. 
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1 Introduction 



A 2- dimensional Almost Riemannian Structure (2-ARS for short) is a generalized Riemannian struc- 
ture that can be defined locally by a pair of smooth vector fields on a 2-dimensional manifold M, 
satisfying the Hormander condition. These vector fields play the role of an orthonormal frame. 

Let us denote by A.{q) the linear span of the two vector fields at a point q. Where A.{q) is 
2-dimensional, the corresponding metric is Riemannian. Where A.{q) is 1-dimensional, the corre- 
sponding Riemannian metric is not well defined, but thanks to the Hormander condition one can 
still define the Carnot-Caratheodory distance between two points, which happens to be finite and 
continuous. 

2-ARSs were introduced in the context of hypoelliptic operators |19| I20| . they appeared in 
problems of population transfer in quantum systems |12| [T3t I14j . and have applications to orbital 
transfer in space mechanics |8l [9]. 2-ARSs are a particular case of rank- varying sub- Riemannian 
structures (see for instance [71 [22l 134]). 

Generically (i.e for an open and dense subset of the set of all 2-ARSs, in a suitable topology), 
the singular set Z, where A.{q) has dimension 1, is a 1-dimensional embedded submanifold and 
there are three types of points: Riemannian points, Grushin points where A.{q) is 1-dimensional and 
dim(A((;) -|- [A,A](g)) = 2 and tangency points where dim{A{q) -\- [A, A](g)) = 1 and the missing 
direction is obtained with one more bracket. One can easily show that at Grushin points A((/) is 
transversal to Z. Generically, at tangency points A.{q) is tangent to Z and tangency points are 
isolated. Normal forms at Riemannian, Grushin and tangency points were established in [3] and 
are described in Figure 1. 

2-ARSs present very interesting phenomena. For instance, geodesies can pass through the sin- 
gular set, with no singularities even if all Riemannian quantities (as for instance the metric, the 
Riemannian area, the curvature) explode while approaching Z. Moreover the presence of a singular 
set permits the conjugate locus to be nonempty even if the Gaussian curvature is always nega- 
tive, where it is defined (see [3|). See also (Sj 131 [TBI IE] for Gauss-Bonnet-type formulas and for a 
classification of 2-ARSs from the point of view of Lipschitz equivalence. 

In this paper we study the Laplace-Beltrami operator in a 2-ARS. The main point is that the 
first order terms explode as a consequence of the explosion of the area on the singular set Z. In 
particular we are interested in its self-adjointness. This is crucial to understand the evolution of a 
quantum particle and of the heat flow in a 2-ARS. 

On M \Z, this operator is defined as the divergence of the gradient. The almost- Riemannian 
gradient can be defined with no difficulty since it involves only the inverse metric which is well defined 
(by continuity) even on the singular set. More precisely, on a Riemannian manifold {A-I,g), V is 
the unique operator from C°°{M) to Vec(M) satisfying 5g(V (j){q),v) = d(j)q{v), V g G M, v G TgM. 
In the two-dimensional case, if in an open set Q C M an orthonormal frame for g is given by two 
vector fields Xi and X2, then we have on 0, 

V(0) =Xi(</.)Xi+X2(</.)X2, (1) 

where by Xi{(j)) we mean the Lie derivative of (j) in the direction of (i = 1, 2). This last formula 
can be used to define the gradient even where Xi and Xi are not linearly independent. In this 
paper, we use formula ^ to define the gradient of a function on a 2-ARS. 

However the divergence of a vector field requires a notion of area and in a 2-ARS the natural 
area (i.e. the Riemannian one) explodes while q approaches the singular set. Let us recall that the 
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Figure 1: The local representations established in |3]. Tangency points are the most difficult to 
handle due to the fact that the asymptotic of the distance is different from the two sides of the 
singular set (see [101 115)). 



divergence of a vector field X on a Riemannian manifold (M, g) is the unique function satisfying 
diwXdijJ = Lxdu where du is the Riemannian volume form. In coor dinates divX = ^ di (^^^') , 
where uj = yjdetg. 

To see how the Riemannian quantities (and in particular doj) explode while approaching the 
singular set Z, choose a local orthonormal frame for the 2-ARS of the form Xi = (1,0) and 
^2 = (0, /(a^, 2/)), where / is some smooth function. This is always possible, see for instance [3]. 
On M \ Z, for the metric g, the area element dio, the curvature K and the gradient of a smooth 
function one has: 

(' ? V 

1 

777 ^dxdy, 

\fix,y)\ 

f{x,y)dlf{x,y)-2dJ{x,yY 



9{x,y) = 

duj{x,y) = 

K{x,y) = 

V(</.) = 
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and for a vector field Y = {Y^, Y'^) one has 

= a,Fi + a,r2-Mri-My2 (2) 

which is not well defined on Z even on vector fields which are linear combinations of Xi and X2. 
In particular it is not well defined on Xi. As a consequence the Laplace Beltrami operator presents 
some singularities in the first order terms. (The choice of the area does not affect the principal 
symbol of the operator). More precisely, one has 

A((^) = d^^ + fd^y^ - ^d,Ct> + f{dyf){dyCP). (3) 

The simplest case is the well known Grushin metric, which is the 2-ARS on for which an 
orthonormal basis is given by, 

X.^d). X.-('>]. (4) 



0/ ' ^ I 

/ \ ^ ^ 

Here the singular set Z is the y axis and on \ Z the Riemannian metric, the Riemannian area 
and the Gaussian curvature are given respectively by: 

9=(\ 1 ), duj = -^dxdy, K = — \. (5) 
It follows that the Laplace Beltrami operator is given by 

:= div(V(</.)) = [dl + x^dl - h^)4>, (6) 

which is the basic model of Laplace Beltrami operator on 2-ARSs and whose self-adjointness will 
be studied in this paper. Notice that the operator "sum of squares", 

A = Xl+Xl = dl + x^dl (7) 

has been deeply studied in the literature. This operator (which in this case is the principal symbol 
of A) can be obtained as the divergence of the gradient by taking as area, the standard Euclidean 
area on R^. However this operator does not give any information about A because of the diverging 
first order term. Moreover in 2-ARS an intrinsic Laplacian with non-diverging terms cannot be 
defined naturally, since an intrinsic area which does not diverge on the singular set is not known. 
Here by intrinsic we mean an area which depends only on the 2-ARS and not on the choice of the 
coordinates and of the orthonormal frame. 

Notice that even if all Riemannian quantities are not defined on M\ Z, classical geodesies do. 
See Section 13.11 for the explicit expressions of geodesies in the Grushin plane. Another interesting 
feature of the Grushin plane, is the fact that a bounded open set intersecting Z has finite diameter 
but infinite area. 

The main purpose of this paper is to study the following question: 
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[Q] Let M be a 2-D manifold endowed with a 2-ARS, and let A be the corresponding Laplace- 
Beltrami defined on C^{M \Z). Is A essentially self-adjoint? 

Notice that a priory one expects a negative answer to this question, since, as explained below, a 
positive answer would imply that neither the heat flow, neither a quantum particle can pass through 
Z, while classical geodesies cross it with no singularities. 

Our main result is an unexpected positive answer to the question [Q] in the case in which there 
are no tangency points and in the case in which M is compact. More precisely we have the following. 

Theorem 1. Let M he a 2-D compact orientahle manifold endowed with a 2-ARS. Assume that 

[HA] the singular set Z is an embedded one- dimensional submanifold of M; 

[HB] for every q G M, A{q) + [A, A]{q) = TgM. 

Let duj be the corresponding Riemannian area and A the corresponding Laplace- Beltrami operator, 
both defined on M\Z. Then we have the following facts 

1. A with domain Cq^{M \ Z) is essentially self-adjoint on L'^{M,dijj). 

2. The domain of A (the self-adjoint extension denoted by the same letter) is given by 

L>(A) = {u G L^{M,duj) \AD,gU G L^{M,duj)] 
where A^^gU is Au seen as a distribution in Q = M \ Z. 

3. The resolvent ( — A + 1)^^ is compact and therefore its spectrum is discrete and consists of 
eigenvalues with finite multiplicity. 

Hypothesis HA is generic (see [3]). Hypothesis HB implies that every point is either a Rieman- 
nian point or a Grushin point (see Figure 1). The hypotheses that M is compact and that there 
no tangency points are technical and the same result should hold in much more general situations. 
The orientability of the manifold can be weakened. Indeed it is only necessary that each connected 
component of Z (which is difl'eomorfic to S^) admit an open tubular neighborhood diffeomorfic to 
(-1,1) X 5^ See Proposition El and El 

Remark 1. Notice that 2-ARSs satisfying HA and HB do exist. See for instance where such 
a 2-ARS has been built on compact orientable manifolds of any genus. Notice moreover that a 2- 
ARS satisfying HA and HB are structurally stable in the sense that small perturbations of the local 
orthonormal frames in the C°° norm, do not destroy conditions HA and HB. 

Theorem [T] has a certain number of implications. First it implies that a quantum particle, local- 
ized for t = on a connected component of M \ ^, remains localized in this connected component 
for any time t G R. The same phenomenon holds for the wave or the heat equation. Indeed, the es- 
sential self-adjointness means that our operator can be naturally and uniquely extended (by taking 
its closure) to a self-adjoint operator without adding any additional boundary condition. But, in 
our case, one possible extension is to take an extension for the operator defined on each connected 
component (for instance the Friedrichs extension) and to "concatenate" them. This possible exten- 
sion is the one which separates the dynamics. For instance this is what we obtain when one defines 
the Friedrichs extension of the usual Laplace operator on R defined on function in (7q°(R \ {0}). 
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We obtain two distinct dynamics with Dirichlet boundary condition at 0. Yet, this is of course not 
the only possible dynamic because the operator defined on function C^(R \ {0}) is not essentially 
self-adjoint. So, the essential self-adjointness means that the unique self-adjoint extension of our 
operator is the one that separates the dynamics on each connected components. That means that 
it is not necessary to add any boundary condition: the explosion of the area naturally acts as a 
barrier which prevents the crossing of the degeneracy zone by the particules. 
More precisely we have the following. 

Corollary 1. With the notations of TheoremUl consider the unique solution u of the Schrddinger 
equation (according to the self-adjoint extension defined in the previous theorem), 

( idtu + Au = , . 

\ u{0) = uoe L'^{M,dio) ^ ' 

with uq supported in a connected component Q of M\Z. Then, u{t) is supported in VI for any t > 0. 
The same holds for the solution of the heat or for the solution of the wave equation. 

Second it is well known that in Riemannian geometry one can relate properties of the Riemannian 
distance to those of the corresponding heat kernels. For instance we have that 

Theorem 2 (Varadhan, Neel, Stroock). Let M he a compact, connected, smooth Riemannian man- 
ifold and d the corresponding Riemannian distance. Let PtiQi,Q2) be the heat kernel of the heat 
equation dt4) = ^Ac/i, where A is the Laplace- Beltrami operator. Define 

Et(qi,q2) ■= -tlogptiqi,q2) 

We have: 

1) Et{qi,q2) — )• ^d{qi,q2)'^, uniformly on M x M as t — )• 0, 

2) Let Cut{qi) C M the cut locus from qi. Then q2 ^Cut{qi) if and only i/ limt^o 52) = 

^V^(i(gi, 92)^, while q2 &Cut{qi) if and only z/ limsupj_j.o ||V^£^t(gi, 52)!! = 00, where \\.\\ is 
the operator norm. 

The first result is due to Varadhan |33| , the second one to Neel and Stroock |28| [29] . 

Both these results do not depend upon the fact that one is using the Laplace-Beltrami operator 
or the principal symbol of the operator in a fixed frame to construct the heat kernel. 

The situation is very different in 2-ARG. Both these results are false for the Laplace Beltrami 
operator. Indeed the distance between two points belonging to two different connected components 
of M\Z is finite while if q belongs to a connected component W of M\Z then ptiq, •) is supported 
in W. 

However a result in the spirit of the one of Varadhan has been obtained by Leandre in |25) 
for the operator "sum of squares". Hence the Laplace-Beltrami operator defined above, has quite 
different properties with respect to the operator sum of squares. The last one is not intrinsic, but 
however keeps tracks of intrinsic quantities as the almost-Riemannian distance. In particular the 
corresponding heat fiow crosses the set Z which is not the case for the Laplace-Beltrami operator. 
For other relations among the heat kernel and the distance in sub- Riemannian geometry see [6]. 

To prove Theorem [T] we start by analyzing the Grushin case ([6]). We first compactify in the y 
variable by considering it on R^; x Ty. 
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By setting / = -y/jxf^ we are reduced to study the essential self adjointness of the following 
operator considered on with the usual euclidian metric 

L = dl + ^'dl-\^,. (9) 

If for a moment we forget the term x^dy which is not relevant for the selfadjointness of this operator 
(and becomes a non-positive potential —k'^x^ after performing Fourier transform in y), we are 
reduced to study the operator 9^ — which is well known in the literature. Indeed we have the 
following result 

Proposition 1. The operator —d'^ + defined on -L^(]0,+oo[) with domain C^(]0,+oo[) is es- 

3 

4- 



sentially self-adjoint if and only c > ^ 



There are different proof of this result in the literature (see Chapter X for an introduction) , 
often leading to some stronger statement (for instance for potential V > ■^)) and generalizations 
to higher dimensions. 

The rest of the proof for an almost Riemannian structure consists in generalizing this result for a 
normal form around a connected component of the singular set. The main tools are Kato inequality 
and perturbation theory. 



Remark 2. Notice that, as remarked in the introduction of ]21^. solutions of the Schrddinger equa- 
tion, with the Laplacian defined as sum of square d"^ + x^dy, display a total lack of dispersion. It 
would he interesting to understand if the same holds for the Laplace- Beltrami operator defined by 
It would be also interesting to study the behaviour of a wave packet moving towards the singular 
set, to understand if there is a reflection or a dispersion. Since in Theorem U\ we prove that the 
spectrum is discrete, it seems natural to expect a reflection. 

The structure of the paper is the following. In section [2] we briefly recall the notion of almost- 
Riemannian structure. An original result is given in Proposition [21 in which we globalize the normal 
form of type {F2) around a connected component of Z in the compact orientable case and under 
the assumption that there are no tangency points. In Section [3] we analyze the Grushin case, both 
for what concerns geodesies and the self- adjointness of the Laplace-Beltrami operator. The main 
result, is proven in Section HI 

As a byproduct of our studies, we obtain that the Laplace-Beltrami operator is not essentially 
self-adjoint in the case in which an orthonormal basis is given by the vectors (1,0), (0, |x|") with 
a £ (0, 1) (see remark |4|). This fact motivates some remarks on metrics of this kind (which do not 
enter in the standard framework of almost-Riemannian geometry, since the vector fields are not 
smooth). For instance we show that in this case, even if the area element explodes, the area of a 
bounded open set intersecting Z is finite (which is false for 2-ARSs). Moreover we show that there 
exist regular curves passing through Z with a velocity not belonging to the span of Xi and X2, 
having finite length. See Section 13.31 

Finally, in the appendix, we discuss the Martinet case to show that this kind of phenomenon 
appears also in sub-Riemannian structure of constant rank but which are not equiregular i.e. in 
which the grow vector depends on the point (see for instance |27|). In this case the role of the 
Laplace-Beltrami operator is played by the intrinsic sub-Riemannian Laplacian defined via the 
Popp measure (see [2l[27]). While aproaching the Martinet surface the Popp volume explodes and 
the first order coefficient of the sub-Riemannian Laplacian does as well. 
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Theorem 3. Consider the sub-Riemannian structure in M = T^,. xRy xT^ for which an orthonormal 

2 

basis is given by Xi = (1,0, X2 = (0,1,0). Then the corresponding intrinsic sub-Riemannian 
Laplacian which is given by Asr = (Xi)^ + (^2)^ — with domain C^{M\{y = 0}) is essentially 
self-adjoint on L'^{M,duj), where doj = -^dxdydz is the Popp volume. 

2 Basic Definitions 

In this section we recall some basic definitions in the framework of 2- ARSs following [SJ . 

Let M be a smooth surface without boundary. Throughout the paper, unless specified, manifolds 
are smooth (i.e., C°°) and without boundary; vector fields and differential forms are smooth. The 
set of smooth vector fields on M is denoted by Vec(M). The circle is denoted by or T depending 
on the context. Cq° denotes the set of smooth functions with compact support. 

2.1 2-Almost-Riemannian Structures 

Definition 1. A 2-dimensional almost- Riemannian structure (2-ARS) is a triple S = {E,f, {■,■)) 
where E is a vector bundle of rank 2 over M and (•, •) is a Euclidean structure on E, that is, (•, ■)q 
is a scalar product on Eq smoothly depending on q. Finally f : E ^ TM is a morphism of vector 
bundles, i.e., (i) the diagram 

E — ^TM 

TT 

M 

commutes, where vr : TM — t- M and tte : E ^ M denote the canonical projections and (ii) f is 
linear on fibers. 

Denoting by T{E) the C°° {M)-module of smooth sections on E, and by f^, : T{E) — )■ Vec(M), 
the map a 1— )■ f*{cr) := fo a. We require that the submodule of Vec{M) given by k = f*(r(£')) to be 
bracket generating, i.e., Lieq{k) = TqM for every q G M. 

Here Lie(A) is the smallest Lie subalgebra of Vec(M) containing ▲ and Lieq(A) is the linear 
subspace of TqA-I whose elements are evaluation at q of elements belonging to Lie(A). The condition 
that ▲ satisfies the Lie bracket generating assumption is known also as the Hormander condition. 

We say that a 2-ARS {E,f, {■,■)) is orientable if E is orientable as a vector bundle. Notice 
that one can build non-orientable 2-ARSs on orientable manifolds and orientable 2-ARSs on non- 
orientable manifolds. See [3] for some examples. We say that a 2-ARS {E,f, {■, •)) is trivializable if 
E is isomorphic to the trivial bundle M x R^. A particular case of 2-ARSs is given by Riemannian 
surfaces. In this case E = TM and f is the identity. 

Let S = {E,f, {■,■)) be a 2-ARS on a surface M. We denote by k{q) the linear subspace 
{V{q) \V € a} = f{Eq) C TqM. The set of points in M such that dim(A(g)) < 2 is called singular 
set and denoted by Z. Since A is bracket generating, the subspace A{q) is nontrivial for every q 
and Z coincides with the set of points q where A is one-dimensional. 

The Euclidean structure on E allows to define a symmetric positive definite C°°(M)-bilinear 
form on the submodule A by 

G : A X A ^ C°°(M) 
G{V,W) = {av,aw) 
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where ay, are the uniquqj sections of E satisfying f o cry = V^, f o aw = W . 

At points q £ M where fl^;^ is an isomorphism, i.e. on M \ G is a tensor and the value 
G{V, W)\q depends only on V{q)^ W{q). In this case G defines a Riemannian metric g via 

gg{v,w) = G{V,W){q), v,weTgM, 

where V and W are two vector fields such that v = V{q) and w = W{q). 
This is no longer true at points q where fl^;,^ is not injective. 

If (o"i,(72) is an orthonormal frame for (•,•) on an open subset Q of M, an orthonormal frame 
for G on is given by (j^ o o"!, f o (72). Orthonormal frames are systems of local generators of ▲. 
For every q £ M and every v £ A{q) define 

Gg{v) = inf{('u,'u)g I u G Eg,f{u) = v}. 

An absolutely continuous curve 7 : [0, T] — )■ M is admissible for S if there exists a measurable 
essentially bounded function 

called control function such that 'y{t) = f{u{t)) for almost every t £ [0,T]. Given an admissible 
curve 7 : [0,T] — )• M, the length of j is 

^(7) = 1^ ^jG^itm)) dt. 

The Carnot-Caratheodory distance (or sub- Riemannian distance) on M associated with S is defined 
as 

d{qo,qi) = inf{^(7) | 7(0) = go,7(^) = Qi^l admissible}. 

It is a standard fact that ^(7) is invariant under reparameterization of the curve 7. Moreover, 
if an admissible curve 7 minimizes the so-called energy functional E{'y) = Jq G^(t)(7(t)) dt with T 

fixed (and fixed initial and final point) then v = ^jGry(^f){j{t)) is constant and 7 is also a minimizer 
of £{■). On the other hand a minimizer 7 of £{■) such that v is constant is a minimizer of E{-) with 

T = e{^)/v. 

The finiteness and the continuity of d{-,-) with respect to the topology of M are guaranteed by 
the Lie bracket generating assumption on the 2-ARS (see [5]). The Carnot-Caratheodory distance 
endows M with the structure of metric space compatible with the topology of M as differential 
manifold. 

When the 2-ARS is trivializable, the problem of finding a curve minimizing the energy between 
two fixed points qo,qi G M is naturally formulated as the distributional optimal control problem 
with quadratic cost and fixed final time 

2 T 2 

q = ^UiXi{q), men, / (t) dt min, q{0) = qQ, q(T) = qi. 

^=l Jo 

where {Xi,X2} is an orthonormal frame. 

^the uniqueness is consequence of the fact that we assume E of rank two and ▲ Lie bracket generating. 
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2.2 Geodesies 



A geodesic for 5 is a curve 7 : [0, T] — )■ M such that for every sufficiently small nontrivial interval 
[^11^2] C [0,T], j\[t-j^fy] is a minimizer of E{-). A geodesic for which G^(j)(7(t)) is (constantly) equal 
to one is said to be parameterized by arclength. The local existence of minimizing geodesies is a 
standard consequence of Filippov Theorem (see for instance [5]). When M is compact any two 
points of M are connected by a minimizing geodesic. 

Locally, in an open set $7, if {Xi,X2} is an orthonormal frame, a curve parameterized by 
arclength is a geodesic if and only if it is the projection on of a solution of the Hamiltonian 
system corresponding to the Hamiltonian 

H{q,X) = ^{iXX,{q))^ + iXX2{q))^), q G Q, X G T*n. (10) 

lying on the level set H = 1/2. This is the Pontryagin Maximum Principle |30| in the case of 2-ARSs. 
Its simple form follows from the absence of abnormal extremals in 2-ARSs, as a consequence of the 
Hormander condition see [3]. Notice that when looking for a geodesic 7 minimizing the energy 
from a submanifold T (possibly of dimension zero), one should add the transversality condition 

A(o)r^(o)r = 0. 

2.3 Generic 2-ARSs and normal forms 

A property (P) defined for 2-ARSs is said to be generic if for every rank-2 vector bundle E over M, 
(P) holds for every f in an open and dense subset of the set of morphisms of vector bundles from E 
to TM, endowed with the C°°- Whitney topology. 

Define Ai = A and A^+i = Afc -|- [A, A^]. We say that S satisfies condition (HO) if the following 
properties hold: (i) Z is an embedded one-dimensional submanifold of M; (ii) the points q £ M at 
which A2('?) is one-dimensional are isolated; (iii) A3 ((7) = TqM for every q S M. It is not difficult 
to prove that property (HO) is generic among 2-ARSs (see [3j). This hypothesis was essential to 
show Gauss-Bonnet type results for ARSs in [3l HI EZ]- The following theorem recalls the local 
normal forms for ARSs satisfying hypothesis (HO) (see Figure 1). 

Theorem 4 ([3]). Consider a 2-ARS satisfies (HO). Then for every point q £ M there exist a 
neighborhood U of q and an orthonormal frame {Xi,X2} of the ARS on U such that, up to a 
change of coordinates, q = (0, 0) and {Xi,X2} has one of the forms 

(Fl) X,ix,y) = i:, X^ix.y) = e^^^'V) 

(F2) Xi(x,y) = ^, X2(x,y) = xe'^(-'J')^, 

(F3) Xi(x,y) = |:, X2{x,y) = {y - x^^{x))e^^-'y) ^, 

where (p, ip and ^ are smooth functions such that 0(0, y) = and tp{0) 7^ 0. 

Let 5 be a 2-ARS satisfying (HO). A point q € M is said to be an ordinary point if A(g) = TqM, 
hence, if S is locally described by (Fl). We call q a Grushin point if A((/) is one-dimensional and 
A.2{q) = TqM, i.e., if the local description (F2) applies. Finally, if A{q) = A2{q) has dimension one 
and ^siq) = TqM then we say that g is a tangency point and S can be described near q by the 
normal form (F3). 

Notice that under hypotheses HA and HB of Theorem [H HO is fulfilled and there are no 
Tangency points. 
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In the compact case, the fohowing proposition permits to extend the normal form (F2) to a 
neighborhood of a connected component of when there are no tangency points. 



Proposition 2. Consider a 2-ARS on a compact orientable manifold satisfying (HO). Let W be a 
connected component of Z containing no tangency points. Then there exists a tubular neighborhood 
U of W and an orthonormal frame {Xi,X2} of the 2-ARS on U such that, up to a change of 
coordinates, W = {(0, y), y G S^} and {Xi,X2} has the form 

Xi{x,y) = ^, X2(x,?/)=xe^(-'^)^. (11) 

Proof. The proof consists in using as first coordinates the distance from W and it is very similar to 
the proof of Theorem HI which contains a local version of Proposition [2] (see [3] pp. 813-814, proof 
of Lemma 1 and proof of Theorem 1 for Grushin points). Here we just explain which modifications 
are necessary. 

Since M is assumed to be compact then W is diffeomorphic to S^. Moreover, since M is 
assumed to be orientable, a sufficiently small open tubular neighborhood of W is diffeomorphic to 
(— 1, 1) X S^. Then consider a smooth regular parametrization 3 a ^ w[a) of W. 

Let a I— 7- Ao(a) G T^^^-^M be a smooth map satisfying H{XQ{a),w{a)) = 1/2, where H is the 
Hamiltonian of the PMP (llOp and Ao(a) -L T^(q)VF. Notice that such a map exists since we are 
assuming that M is orientable. 

Let E{t, a) be the solution at time t of the Hamiltonian system given by the Pontryagin Maxi- 
mum Principle with initial condition {q{0), \{0)) = (t(;(a), Ao(a)). With the same arguments given 
in the proof of Lemma 1 in j3], one shows that E(t, a) is a local diffeomorphism around every point 
of the type (0, a), q G S^. Using the fact that E is a global diffeomorphism from {0} x toW and 
by suitably reducing e one gets that it is a diffeomorphism between (— e, e) x (for some e > 0) 
and a tubular neighborhood of W. This permits to use as coordinates in U the pair (t,a). Indeed 
t is the distance from W which we proved to be smooth in U. As in [3], one builds an orthonormal 
frame in U defining the vector field Xi by 

Xi{t,a) =dtE{t,a) =dt. 
As in [3] one obtains that the second vector field of the orthonormal frame has to be of the form 

Remark 3. Notice that in Proposition^^ the hypothesis that M is orientable can be weakened by 
requiring that W admits an open tubular neighborhood diffeomorphic to (—1, 1) x 

3 The Grushin case 

As already mentioned the Grushin plane is the trivializable almost Riemannian metric on the (x, y) 
plane for which an orthonormal basis is given by 

In the sense of Definition [T] it can be seen as a triple {E, f, {-,■)) where E = xR^, f((x, y), (o, b)) = 
{{x,y), {a,bx)) and (•, •) is the standard Euclidean metric. 
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Figure 2: Geodesies and front for the Grushin plane, starting from the singular set. 



3.1 Geodesies of the Grushin plane 

In this section we recall how to compute the geodesies for the Grushin plane, with the purpose of 
stressing that they can cross the singular set with no singularities. 

Setting q = {x,y) and A = (Ai, A2), the Hamiltonian (|lUp is given by 

H{x,y,Xi,X2) = ^iXl + x'Xl) (13) 
and the corresponding Hamiltonian equations are: 

X = Ai, Ai = — xA| 

y = x^X2, A2 = (14) 

Geodesies parameterized by arclength are projections on the (x, y) plane of solutions of these 
equations, lying on the level set H = 1/2. We study the geodesies starting from i) a Grushin point, 
e.g. (0,0) ii) an ordinary point, e.g. ( — 1,0). 

Case (x(0),y(0)) = (0,0) 

In this case the condition i?(a;(0), y(0), Ai (0), A2(0)) = 1/2 implies that we have two families of 

geodesies corresponding respectively to Ai(0) = ±1, A2(0) =: a E R. Their expression can be 
easily obtained and it is given by: 



r x{t) = ±t, y{t) = ifa = 

Some geodesies are plotted in Figure 2 together with the "front" i.e. the end point of all geodesies 
at time t = 1. Notice that geodesies start horizontally. The particular form of the front shows the 
presence of a conjugate locus accumulating to the origin. 
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Case {x{0),y{0)) = (-1,0) 

In this case the condition iJ(x(0), y(0), Ai (0), A2(0)) = 1/2 becomes + A| = 1 and it is convenient 
to set Ai = cos(), A2 = sin(), G S^. The expression of the geodesies is given by: 

x{t) = t-l, y{t) = if =0 

xlt) = -t-1, y{t) =0 if = vr 

sinf - 1 smO) 2t - 2 cos() + ^^^i^-^^^f^O) (16) 

sin() 4sin(j 

Some geodesies are plotted in Figure 3 together with the "front" at time t = 4.8. Notice that 
geodesies pass horizontally through with no singularities. The particular form of the front shows 
the presence of a conjugate locus. Geodesies can have conjugate times only after intersecting Z. 
Before it is impossible since they are Riemannian and the curvature is negative. 




Figure 3: Geodesies and front for the Grushin plane, starting from a Riemannian point 



3.2 The Laplace-Beltrami operator on the Grushin plane 

In this section we illustrate in details the results of the paper in the Grushin case. For the study of 
the Laplace Beltrami operator it will be convenient to compactify the y direction considering x G R 
and y G T. 

As mentioned in the introduction, the Laplace-Beltrami operator on the Grushin plane can be 
written as 

We expect that the Laplacian is the Friedrichs extension associated to the positive quadratic form 

Q{u,u) = / |Vnp = / - — -\dxu\'^ + \x\ \dyu\'^ dxdy 
Jr^ Jr^ \x\ 
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Let us make the change of variable / = A/fJcfs' which is unitary from L^(R x T,du}) to L^(R x 

T,dxdy), so that {fi, f2)L^-RxT,d<^) = {9i,92)L^n.xT,dxdy)- 
We compute the operator in the new variable: 



A/ = dlf + x^d^f-j^d^f 



d,g + x dyg- 



x\Lg. 



Hence we are left to study the operator Lg = d'^g + x'^dyg — -^g on L^(R x T, dx dy). Decomposing 
a function g[x^ y) in the Fourier basis in the variable y, we get the decomposition 

L2(RxT)«efgziffc 

where Hj. = I/^(R) and the operator L acts on each Hj. by 

3 

{Lu)k{x) = dluk - x^k^Uk - -^^k = dluk - VkUk- 

with Vk{x) > ^p-. But, we know that in dimension 1, the operator —9^ + V with domain 
C^(]0,+oo[) is essentially self-adjoint on L2(]0,+oo[) ii V > ^ (see [3l], Theorem X.IO for 
the proof of the limit point case at and Theorem X.8 at +00). Hence, each operator d'^Uk — Vk is 
essentially self-adjoint. As a consequence L is essentially self-adjoint as well. 

Remark 4. V = is exactly the limit singular potential to have —d'^ -\-V in the limit point case 
at zero and hence essentially self-adjoint (see ]31\ Theorem X.IO, page 159]). It is surprising to find 
the same constant | for the Grushin operator. 

Notice that if we consider as orthonormal basis of the metric the vector fields X = and 
y = ^1 ^]^a^ with a > 0, we obtain for the Laplace- Beltrami operator, 

Au = d'iu -I- |xP°(9^u — —dxU. 

^ X 

Making the change of variable u = |x|°/^f (which is unitary from L^(Rx T, dio) to L^(Rx T, dx dy) ) 
we get for the transformed operator 

~^v = dlv-l[l + l)^^+\x\--dlv, 

which is essentially self-adjoint if and only if a £ [l,-|-oo[. 



3.3 'VGnishin 

In this section, motivated by Remark HI we briefly discuss the generalized Riemannian structures 
on R^ for which an orthonormal basis is given by 

^i=( J), ^2=( 1° ), aG(0,l). (17) 
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Notice that this generahzed-Riemannian structure does not fit the definition of almost- Riemannian 
structures, since the vector fields are not smooth. As usual call Z the set in which the two vector 
fields are not linearly independent. It is interesting to notice the following facts, which in some 
cases are very different with respect to the standard Grushin case: 

• The Lie bracket between X\ and is not well defined on Z. Indeed [Xi, X2] = (0, a|x|^"~"^^)"^. 
However it is not difficult to prove that the corresponding distance is continuous and endows 

with its standard topology (note that it is not the case if a > 0). 

• The curvature is given hy K = — and the area element by duj = -^dxdy. Hence the 
area element explodes, but the area of a bounded open set intersecting Z is finite (and its 
diameter is finite as well). This is certainly the reason why the Laplace-Beltrami operator 
is not essentially self-adjoint and, as a consequence, a quantum particle or the heat flow can 
pass through the set Z. Notice that for the Grushin plane a bounded open set intersecting Z 
has finite diameter but infinite area. 

• There exists regular curves passing through Z with a velocity not belonging to the span of 
Xi and X2 and having finite length. For instance take a = 1/2 and the curve 7(t) := (t,t), 
t £ [-1,1]. We have 7(0) ^ span{Xi(7(0)),X2(7(0)} but its length is fini te. Inde ed it 
corresponds to the controls ui = 1 and U2 = l/|t|^/^. Hence ^(7) = J^-^ \J\ -|- = 
2(-\/2 -|- arcsinh(l)) ~ 4.59. Such phenomenon does not appear on the singular set of a 2-ARS, 
as a consequence of the smoothness of orthonormal frames. 

4 Laplace-Beltrami on an Almost-Riemannian structure 

4.1 Computation in local coordinates 

The aim of this subsection is the proof of the following theorem. 

Theorem 5. Let g a Riemannian metric onQ = {x = {x,y) £ H x T \ x ^ 0} of the form diag(^l, x~'^e 
where ^ is a smooth function which is constant for large \x\. Let A be the corresponding Laplace- 
Beltrami operator with domain Cq°(0). 

Then, A is essentially self-adjoint for L^(R x T, dhj), where doj is the corresponding Riemannian 
area on Q.. 

The same result holds by replacing Q with Q± where r2_ = {x = {x,y) £ H x T \x < 0} and 
n+ = {x = {x,y) G R X T|x > 0}. 

A simple consequence of this theorem is that the only self-adjoint Laplace operator that can be 
constructed by extension of the natural one (i.e. the Laplace-Beltrami operator defined on C^(J7)) 
preserves the decomposition L'^{Q,,duj) = L'^[U^,duj) ©-*- L'^[Q.^,du) . 

Indeed, if we denote A~ (resp A^) the unique self-adjoint extension of the natural Laplace 
operator with domain Cq°{Q-) (resp CQ°(r2+)), we can define the operator A{u- +u+) = A~u_ -|- 
A^n+ with the natural domain inherited from the one of A~ and A+. The operator A is a self- 
adjoint extension of the Laplace-Beltrami operator with domain CQ°(f2) and it is the unique one by 
Theorem [5l Another way of seeing this fact is from the point of view of the evolution equation. 
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Corollary 2. With the notations of Theorem\^ consider the unique solution u of the Schrddinger 
equation (according to the self-adjoint extension defined in the previous theorem), 



idtu + Au = 

n(0) = uo G L2(R X T,dw) 



(18) 



with Uq supported in 0+. Then, u{t) is supported in 17+ for any t > 0. The same holds for the 
solution of the heat or for the solution of the wave equation. 

This corollary is a simple consequence of the essential self-adjointness of the Laplace-Beltrami 
operator as discussed above. 

This is in strong contrast with the classical dynamics associated with this metric where the 
geodesies can "cross" the barrier {x = 0} (see Figure 3). The proof of Theorem [5] relies on a change 
of variable similar to the one we did for the Grushin operator which leads to an operator that can 
be written as "Grushin type + singular potential". Then, we follow the proof of Kalf and Walter |23) 
(which was himself inspired by B. Simon [32]) using Kato inequality for the main part. The other 
part can be treated by perturbation theory. Note that some related results for singular potential 
were proven in |18) and |26| . 



Proof of Theorem\^ In coordinates, the Laplace-Beltrami operator has the form 

^ = Wd-.[ug^^d,] 



where we have denoted uj = ^Jdetg the Riemannian volume. 

For u G L^(R x T,duj), let us make the change of variable u = lo~^I'^v which is unitary from 
L^(R X T, dw) to L^(R x T, dx dy) and let us compute its action in the new variable. 

^1/2^^-1/2^ := Av 



UJ 



T7^E[' 



g'^uj{djUJ-^/'^) V + g'^uj{djUJ-^/^)diV + {dig'^)u;^^^djV + g'^ [diU^'^)djV + g'^ {oj^'^)dijV 
'g'^u{djU-^/^)\ V + g'^u{d,u-^''^)djv + {dig'^)u^''^djV + g'\diOJ^/^)djV + g'\uj^/^)dijV 



where we have used the symmetry of g. 
We have: oj{diUJ-'^/'^) + {diuj^/'^) = - 



1 



UJ 



1/2 



di 



g'^ujid.uj-^'^] 



^UJUJ 



E 



'^I'^diUJ + \uj ^I'^diUJ = 0. Moreover let 



2a;i/2 



ci\d,uj)uj-^l^) 



E 



-^a,(5*^)9,(lna;) - ^9^' + -^J\d,uj){djUj) 
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Hence we get 



Av = ^g'^dijv + {dig'^)djv + ^xv 



= diVeCVv) + fIV 

= diVe{G~^V eV) + fIV . 

In the last formula we have used matrix notation where = {g^^) and divg, Vg are the divergence 
and the gradient of the Euclidian space. Now let us specify the computation in our diagonal metric. 



ijj 
dioj 



1 



sgn(^^_$ 



-r-re , 



X 



(52$) 



e , 



^2$ 

— i — r' 

|x| \x\ 

j_ 2sgn(x)ai$ (9i$)2 _ 

I I ^ 119 ^ 11^ 11^' 



1 sgn(2;)ai$ (9i$)^ Sf^* 



+ 



1 sgn(x)c?i<I> ^ ((?i<I>)2 



4|xP ■ 2|x| ' 4 



I 2 2$ 

+x e 



4 X 



I 2 2<E> 

+x e 



sgn(x)gi$ _ {di<^>f dj^ 

2\^\ 4 ^ ~ir 



41x1 



Now, we are able to conclude by the following two Lemmas. The first one proves that the main 
part without ^2 is essentially self-adjoint while the second one treats ^2 as a perturbation. 

Lemma 1. Let G given by (i^ag(l,x2e2*(^'2')) with |$| <C forG >0. 

The operator L : v ^ — diVe{G~^'V ev) + 4|fj2^ defined on L'^(R x T, dx dy) with domain Cq°{Q) 
is positive and essentially self-adjoint. 



Hence, since L defines a positive symmetric operator, its unique self-adjoint extension is the 
Friedrichs extension (in the following still denoted by L). 

Lemma 2. The operator of multiplication by ^2 (defined with domain G^{Q))) is infinitesimally 
small with respect to L. So, by the theorem of Kato-Rellich (Theorem X.12 of ]31^). L + ^2 is 
essentially self-adjoint. 

To conclude the proof of Theorem [5] we are left to prove the two lemmas. 
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Proof of Lemma [IJ We follow closely the proof of Kalf and Walter |23| . Actually, we mimick the 
proof of the 1-D case and we use the fact that our operator is the classical Laplacian with potential 
for function only depending on x: —df + It is enough to prove that Ran(L + /) is dense in 

L^, see Theorem X.26 of |31) . So, let h £ L^(0) such that h _L Ran(L + /) and let us prove h = 
(we can assume h real valued without loss of generality). 



Set 



l+|a;|3/2 



We have (3{x) > and /3 G L^{n). We obtain, for x / 0, 



-diVe(G-iVe/5) + 



4 b 



r/3 + /3 



4(l + |x|3/2)3 



/3 + /3 



Let r/, C G C°°([0,+oo[) so that 



C(s) 



if s < 1/2 

1 if s > 1 
1 if s < 1 
if s > 2 



Denote /3„(x) = r/(n|3;|)C(^)/3(x) G C^{^) so that L/3„ is well defined in L'^^Q). Moreover, we 
have 



/3(x) = 0(e-N), 



|V/3(x)| = as x ^ 

|V/3(x)| =0(e~l^l), asx^oo 



so, L/3„ + /3n = ??(n|x|)C(^)X - 25i f?(n|x|)C(^) dif3 - df r/(n|x|)C(^) /? converges weakly in 
to X G L^(0) such that X{x) > 0. So, we can write 

{X,\h\)^, = hm (L/3„ + /3„,|/i|)^2 >0 

n— >oo 

But, we will use the Kato inequality: 
Lemma 3 (Kato's inequality, see Lemma A of |24|). 

For V G Lj^^{Q) real valued such that diVe{G~^V ev) G Lj^^{Q), we have, in the sense of distributions 
on 

-diVe{G-^Ve\v\) < -{sgnv)diVe{G~^Vev). 



Notice that the Kato inequality of |24| can only be applied on il^ := {(x, ?/) G O : |x| > e} 
where the metric is positive definite. This does not create difficulties since we are applying this 
Lemma in the sense of distributions on il.. 

Applying the Kato Lemma to h and /3„ > as test function, we get 



So iX,\h\) 



L2 



{Lf3n + /3„, \h\)^2 < ifSn, {sgnh){Lh + h))^2 < 0. 
and h = 0. 



□ 
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Proof of Lemma The quadratic form associated to L + 1 is 



Q{u,u) = [ (G ^Vgti • Vgu) H — r^-TT InP + ItiP dxdy 

Jn M^r 

For u G C^(r2), and by Cauchy-Schwarz inequality, we have 

Q{u,u) = {Lu,u)j^2 + 11^11^2 < ||-^^w||^2 \\u\\]^2 + ||'u||^2 < e ||i^^||i2 + + 1^ 
Moreover, since $ is smooth and constant for large x, we have 

Ii/i2ii|ii2 < c 



2 

L2 



< CQ{u,u) <Ce||Lu||^2 + (C^ + l)||n||^2. 



This yields the result since e is arbitrary small and the estimate can be extended to any u G D{L) 
by taking closure. □ 

The proof of Theorem [5] is concluded. □ 

The next Lemma is the first step to establish the compactness of the resolvent of A in the 
compact case. Of course, since here we are dealing with the case of R x T, we can not expect 
compactness without adding a cut-off function. 

Lemma 4. Denote by A the positive self-adjoint operator defined by Theorem\^ Then, the truncated 
resolvent p{x){—A + where p G Cg°(R), is well-defined and compact from L^(R x T,duj) to 

itself. 

Proof. The fact that (—A + is well defined comes from the positivity of —A. For the com- 
pactness it is equivalent to prove that the operator p{x){—A + is compact when defined on 
L^(R X T,dxdy) with the metric induced by the Lebesgue measure. 

We begin by showing that it is the case for L defined above. Let Un be a real valued sequence 
in L^(R X T,dxdy) with norm bounded by 1. Denote := p{x){L -\- V)~^Un ■= p{x)fn. In 
particular, Q{fn, fn) = {{L + l)/„, /„) = (tin, fn) = Jq G-'^Vefn ■ ^efn + ^ \ fn? + |/nP dx dy is 
bounded in L^(R x T.dxdy). Set e > 0. Take ip £ C^(R), supported in ]- e,e[ with < ip < I 
and 99 = 1 in a neighborhood of 0. Split /„ = ip{x)fn -|- (1 — (p{x))fn ■= fn,i + fn,2- Then, we have 
|x(^(x)|2 < £2 and MnWh = I '"''Si'^" < ^^QifnJn) < Ce\ This gives also < Ce. 

Now that (f is fixed, /„^2 = (1 ~ fix))fn = (1 ~ fix))p{x)fn is supported in a fixed compact of 
and bounded in . Indeed, 

||Ve/„,2|li2 < 2(11(1 - ip{x))p{x)VJn\\h + WM^i^) pi^)) fn\\l2) < CQ{fn, fn) < C 

where the second inequality comes from the fact that there exist a constant C such that Id < CG~^ 
(in the sense of quadratic forms) on the support of (1 — (p{x))p{x). We conclude by invoking the 
compact embedding of into Lp' on compact sets so that /„^2 is convergent up to extraction. 
Actually, we have proved that for any e > 0, we can find an extraction 7 such that /^(n) can be 
written /.^(„) = /^(„),i + f^{n),2 with ||^2 < e and /^(„),2 convergent. By choosing e = |, 
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p G N, and by a diagonal extraction argument, we easily get that we can find a subsequence such 
that /■y(^n) is a Cauchy sequence and converges. This gives that p{L + l)""*^ is compact. It remains 
to prove the same result for L + ^2- 

Again, let be a sequence in L2(R X T,dxdy) with norm bounded by 1. Denote fn '■= 
p{x){L + /i2 + ^)~^'>J'n ■= p{x)fn- Thanks to Lemma[2l we get 

||(L + < \\{L + p2 + l)/n|L2 + \\p2fn\\L^ < + S \\{L + + 

so that we get by absorption for e small enough that (L + !)/„ is bounded. Since we have proved 
that p{L + is a compact operator, we get that /„, = p(x)/„ = p{L + (L + !)/„ is relativelly 
compact. □ 



4.2 Case of a compact manifold 

In this section, using the previous results, we prove the three statements of Theorem [TJ For the 
first statement, we have to prove the following. 

Proposition 3. Consider a 2-ARS on a compact manifold M satisfying hypotheses HA and HB. 
Denote by A the Laplace operator defined for functions of C^{0,) where = M \ Z. Then, A is 
essentially self-adjoint on L'^(M,duj). 

Proof. Let M = Ujg/ ilj be a finite covering of M such that for every connected component W oi Z 
there exists i G / such that W and H fij =0 for j € I , j ^ i and an orthonormal frame for 
the 2-ARS in f]j is given by ([TT]) . Moreover assume that if Qj (j G /) does not contain any Grushin 
point, then an orthonormal frame for the 2-ARS in is given by the normal form (Fl). This is 
possible thanks to Theorem U and Proposition [2] 

Theorem [5] yields the result in local coordinates around a connected component of the singular 
set. We only have to extend it to the whole manifold. 

So, let (^j)jg/ be a partition of unity associated to that is \E', G C^{Qi) and X^ie/ = 1- 

We can also assume that >I'j = or = 1 in a neighborhood of Z, so that V^'j and A^j are C°° 
compactly supported in Q. Here we use the fact that we have global coordinates around every 
connected component of Z. 

So again, let u G L^(il) such that u _L Ran(— A -|- /) and let us prove u = (we can assume u 
real valued without loss of generality). We denote Uj = ^'jU G L^(r2j). 

Since the domain of definition of A contains C^{Q), we get that u is solution of —Au + u = in 
the sense of distributions on Q. By elliptic regularity, we get that u G C°°{0,). So, the only problem 
is the possibility to compute integration by part around the degeneracy points. 

Let rii with i G / be an open set around which we can find a coordinate system so that the 
metric g takes the diagonal form (l, x~^e~^*(^'^)) that we extend arbitrarily on R x T with $ 
smooth constant for large x. We will denote Aj, Vj and diOi the Laplacian, gradient and area 
corresponding to this extension on R x T. Since Ui G L^(r2j), we can consider the function Ui 
in this local coordinates (in what follows we will not distinguish Ui with its representent in local 
coordinates) and make the computation (in the sense of distributions of and of R* x T in local 
coordinates) : 

-AiUi + Ui = -2ViU ■ Vi^i - {Ai^i)u. 
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Remark that in these local coordinates u has only a true meaning for small x but we can then extend 
this equality on R* x T since is compactly supported. Moreover, Vj^'j and Aj'!',, is supported 
outside of the zone of degeneracy Z where u is C°° . So, we get — AjUj + tij G Cq°(R* x T) and so 
— AjUj + tij G L?'{M,doJi). This is true in the sense of distributions, but that means that for any 
if G C^(R* X T), we have 

|(ni,-Ai(/? + 99)i2(R.xT,da;i)l ^ II '/'II L2 (R* xT.do;,) • 

In particular, that means that Uj belongs to the domain of the adjoint of the operator — Aj + I 
with domain C^(R* x T). But Theorem [5] gives that this operator is essentially self-adjoint so, 
Ui G D{Li) where Lj is the self-adjoint (Friedrichs) extension of — Aj -|- / on R* x T. In particular, 
we can write (the right hand side has to be understood as a limit for a sequence Ui^n in C'o°(^) 
converging strongly to Ui for the norm of the quadratic form) 

(Uj, -Aj-Uj nj)i2(R.xT,dwO = l|Vj^j|lL2(R.xT,dwO + ll^«llL2(R*xT,daJi) • 

This quantity also makes sense when Ui is considered on the manifold M and the extension of the 
metric was chosen so that — AjUj + Ui = —Aui + Ui in the sense of ditributions and VjUj = Vuj. 
The same result holds for i £ I corresponding to a Riemannian zone. 

Moreover, if i ^ j € I, the common support of Ui and Uj do not intersect the degeneracy zone 
Z. So, these functions are C°° in this zone and we can write 

{ui, -Auj + %)l2(m) = (Vuj, Vuj)l2(a/) + {ui, Uj)l2(m) 
By summing up, we get 

= (U, -An + u)L2(M,dw) = (Vw, yu)L2(_M,d<^) + {u, u)L2^M,dw) 

and u = 0. □ 

From now on, we will use the same notation A for the self-adjoint extension of the symmetric 
operator A. 

To prove the second statement of Theorem [1] we have to prove the following. 
Proposition 4. The domain of A is given by 

D{A) = {ue L'^{M,dio) \AD,gU G L'^{M,dio)} (19) 

where Ao^gU is Au seen as a distribution inVt = M\Z. 

Proof. Let T be the Laplace operator defined with domain given by (jl9p . We easily see that 
D{T) = {u£ L'^iM,duj)\3C : ^^p e C^{n),\{u, Aip)L2 \ <C\\ip\\^2} = D{A*) = D{A). 

□ 

To prove the third statement of Theorem [T] we have to prove the following. 

Proposition 5. Denote by A the positive self-adjoint operator defined by Proposition^^ then the 
resolvent (—A -|- 1)"-*^ is well defined and compact from L?'{M, duj) to itself with the measure defined 
by the metric. 

Therefore, the spectrum of A is discrete and consists of eigenvalues of finite multiplicity. 
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Proof. The fact that (— A + 1) ^ is well defined comes from the positivity of —A. Now let Un 
be a bounded sequence in L^{M,dLo) and /„ = (—A + l)~^Un- By density, we can assume Un in 
C^{M\Z) and /„ G C°°{M\Z) by elliptic regularity. Then, we have /^^ |V/„|2 + |/„|2 = /„,)^2 
bounded. 

Consider the partition of unity introduced in the proof of the previous Theorem and 

denote := "^ifn- 

Let i S / be an index corresponding to a "Grushin zone". By the formula 

(-A + = ^f.Un - 2V^, • Vfn - iA^^)fn, 

we have that (—A + I)fn,i is bounded in L?'{M). By using Lemma |4l we get that Aj + 

A + I)fn,i is compact for a function p G C^(R) defined in some local coordinate charts. 
Here, we have denoted Aj the Laplacian for an extension of the local metric to R x T. To finish, 
we only have to notice that (— Aj + A + I)fn,i = (— Aj + Aj + l)fn,i = fn,i because 

of the support of (note that it is not the case for (— Aj + because the resolvent depends 
on the extension). 

So, choosing p such that p{x)fn,i = fn,i in local coordinates, we have proved that each sequence 
fn,i is compact for a Grushin zone. The same result holds for a Riemannian zone because the 
equivalent of Lemma H] still holds for any Riemannian extension. By summing up, we get that the 
sequence /„ is convergent, up to extraction, which yields the result. □ 



A The Martinet case 

In this Section, we prove Theorem [3j A sub- Riemannian structure is a triple (M, A,g), where M 
is a smooth manifold, A is a smooth vector distribution of constant rank satisfying the Hormander 
condition, and g is a Riemannian metric on ▲. Let Ai = A and Afc_|_i = A^ + [A, A^]. A sub- 
Riemannian structure is said to be equiregular if the dimension of A^ does not depend on the 
point. 

Almost-Riemannian and sub-Riemannian structures can be treated in the unified setting of 
rank- varying sub-Riemannian structures, see [l] Chapter 3] and [H Definition 2]. Here for sake of 
readability we omit this point of view. 

In this section, we briefly treat the Martinet sub-Riemannian manifold (see for instance |1HI27|). 
defined by M = R^, k{q) = span{Xi{q), X2{q)}, and g{Xi,Xj) = 5ij, i,j = 1,2, where 









■ 







Xi = , X2 = 1 . (20) 





Notice that X3 = [Xi,X2] = | | . Hence Xi, X2 and X3 span the tangent space at any point 

-y 

outside the so called Martinet plane Z = {y = 0}, which is the region where the structure is not 



equiregular. However X4 = [[Xi, X2], X2] = [X3, X2] = I I > hence the Hormander condition is 



1 



fulfilled on the whole space. 
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With this structure, \ 2^ is a 3D sub-Riemannian contact manifold. On such a structure 
it is possible to define intrinsically a volume form which is given by dXi A dX2 A dX^, where 
{dXi,dX2,dX-i} is the dual basis to {Xi,X2,X3}. Such a volume form is independent on the 
choice of the orthonormal basis which define the sub-Riemannian structure, and it is called the Popp 
measure (see |27j and Proposition 8 of [2]). One gets dXi = dx, dX2 = dy and dX^ = ^dx — ^dz 

so that dXi A dX2 A dX^ = —^dx A dy A dz and the corresponding density is doj = ^dxdydz. This 
allows to define the sub-Riemannian Laplacian as the divergence of the sub-Riemannian gradient 
defined as in formula ([T]) (see Remark 14 of [2]), 

A,, = (Xi)2 + {X2f - -X2 = (9. + + ^2 - -dy. (21) 

y 2 ^ y 

Notice that the singularity in the first order term of A^^ appears similarly as it does in the 
Grushin case. We do the same reasoning. 

For simplicity, we compactify in x and z and consider the same structure on M = x Ry x T^. 
We will prove that A^^ is essentially self-adjoint on L^(M, dcj) with domain Cq°{M \ Z). Again, 
we make the change of variable / = y^Jy\g so that 

A../ = VW\ ((5,. + y^d.fg + dig - = ^M^srg. (22) 

So, we are left to prove that A^^ = (dx + \Qz)^ ^^'y~ ip' ^® essentially self-adjoint on L^(T x R x 
T,dxdydz) with domain Cq°{M \ Z). We compute the Fourier transform in x and z and get the 
decomposition 

L\T X R X T) « ©t,z)6Z2%,0 
where Hq^^i-^ = L^(R) and the operator L acts on each H^^i^ i-^ by 

y^ 3 

{Lu)(^k,i){y) = dlu(^k,i) - {k + —lfu(^k,i) - -^^{k,i) = dlu(^k,i) - V{k,i)U(k,i)- 

with V(^k.l){y) — We conclude as in the Grushin case. 

This result suggests the general conjecture that for a sub-Riemannian structure which is rank- 
varying or not equiregular on an hypersurface, the singular set acts as a barrier for the heat flow 
and for a quantum particle. 
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